arXivil504.05627v2 [cond-mat.mes-hall] 23 Jun 2015 


Phosphorene analogues: isoelectronic two-dimensional group-IV monochalcogenides 

with orthorhombic structure 
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The group-IV monochalcogenides SnS, SnSe, GeS, and GeSe form a family within the wider group 
of semiconductor ‘phosphorene analogues’. Here, we used first principles calculations to investigate 
systematically their structural, electronic and optical properties, analyzing the changes associated 
with the reduction of dimensionality, from bulk to monolayer or bilayer form. We show that all 
those binary phosphorene analogues are semiconducting, with bandgap energies covering part of the 
infra-red and visible range, and in most cases higher than phosphorene. Further, we found that they 
have multiple valleys in the valence and conduction band, the latter with spin-orbit splitting of the 
order of 19-86 meV. 


I. INTRODUCTION 

Two-dimensional (2D) materials have been extensively 
studied ever since a monolayer graphene was isolated 
by mechanical exfoliation [I|. Thereafter, the interest 
was promptly extended to other 2D materials, such as 
hexagonal boron nitride (h-BN), layered metal dichalco- 
genides (LMDCs) and phosphorene, to name a few 
yj. In special, phosphorene, a monolayer of black phos¬ 
phorus, adopts an orthorhombic structure different from 
graphene and transition metal dichalcogenides.[l,@ This 
anisotropic structure is in the origin of some of phos- 
phorene’s interesting properties, such as superior flex¬ 
ibility under tensile strain, Q and giant thermoelectric 
coefficient. 0 

However, the waved structure of phosphorene is shared 
by yet another class of 2D materials that has so far eluded 
attention. In the bulk form, group-IV monochalcogenides 
GeS, GeSe, SnS and SnSe all assume structures that can 
be considered derivatives of the orthorhombic black phos¬ 
phorus belonging to the space group Pcmn-D^® (lower 
than black phosphorus, which has only one element and 
therefore belongs to Bmab-D^®). This is the a phase of 
SnS, also known as Herzenbergite, a naturally occurring 
but rare mineral0 

Gurrently the most important prospective application 
of a-SnS is as absorber material for film photovoltaic 
(PV) cells. Although other chalcogenide materials as 
CdTe and GuInGaSe 2 also show high PV efficiencies [lll- 
[T^ . many factors make their usag e difficult such as the 
high cost and toxicity of Gd fTTTIT^. [TB| . In contrast, 
SnS is made of abundant and nontoxic elements, and 
its optical band gap of ^ 1.3 eV[ll|) Si is right in the 
range of the optimal values for solar cells (l.I to 1.5 eV). 
Moreover, according to a recent study, solar conversion 


^ Some of these compounds (at least SnS and SnSe) have a more 
symmetric /3-phase, with space symmetry Cmcm-Dj^, at higher 
temperature.)^ However, according to our calculations, all of 
the four compounds are most stable in the o-phase on the mono- 
layer form. 


efficiencies achieved so far for SnS can be well beyond the 
potential limit for the material due to the poor choice of 
band alignment in the devices. [13 In addition, group- 
IV monochalcogenides may show to be superior to other 
2D semiconductors in properties where anisotropy plays 
an important role, as patent in the record thermoelec¬ 
tric coefficient recently reported for SnSe (ZT=2.6 at 
923 K).[I|| 

Besides, these binary ‘phosphorene analogues’ are ex¬ 
pected to reveal distinct intrinsic properties in monolayer 
form, as some of the lattice symmetry operations, in¬ 
cluding inversion, are only present in bulk and in even- 
numbered layer systems. In this sense they are differ¬ 
ent from phosphorene , [13 where inversion symmetry pre¬ 
vents spin-orbit (SO) splitting. In contrast, as shown in 
the present article, group-IV monochalcogenide monolay¬ 
ers show a large intrinsic spin-orbit splitting at valence 
and conduction band valleys. However, even though a 
few theoretical and experimental works have reported on 
the electronic and optical properties of monolayer or few 
layer SnS,[ll|j US] monolayer properties of this group re¬ 
main poorly explored. 

In this work, we use first-principles calculations to 
investigate electronic, structural and optical proper¬ 
ties of the four aforementioned group-IV monochalco¬ 
genides MX, with M=(Sn, Ge) and X=(S, Se), in the 
phosphorene-like a-phase. We compare the properties 
of monolayer and bilayer with those of bulk for each of 
these materials, highlighting the differences in the elec¬ 
tronic and optical properties. 


II. METHODS 

We use first-principles calculations based on density- 
functional theory to obtain the electronic, structural and 
optical properties of monochalcogenides. We employ a 
first-principles approach based on Kohn-Sham density 
functional theory (KS-DFT)[2j, as implemented in the 
Quantum ESPRESSO code.|22||. The exchange corre¬ 
lation energy is described by the generalized gradient ap¬ 
proximation (GGA) using the PBE[2^ functional. Inter- 



2 


actions between valence and core electrons are described 
by Troullier-Martins pseudopotentials [l^. The Kohn- 
Sham orbitals were expanded in a plane-wave basis with 
a cutoff energy of 70 Ry, and for the charge density, a 
cutoff of 280 Ry was used. The Brillouin-zone (BZ) was 
sampled using a T-centered 1 x 10 x 10 grid following the 
scheme proposed by Monkhorst-Pack[^ . For the opti¬ 
cal properties (dielectric constant and conductivity), a 
finer 1x40x40 grid is employed. The calculation of the 
spin-orbit splitting was performed using noncolinear cal¬ 
culations with fully relativistic pseudopotentials. 

In addition, a hybrid functional approximation for the 
exchange-correlation term (HSE06) [26l| is employed in 
order to give reliable results for the gap energies, which 
are well known to be underestimated when employing 
semilocal GGA approximations. For the hybrid func¬ 
tional bandstructure calculations, we used the Vienna 
ab initio simulation package (VASP) [2^ with the 
projector-augmented wave potentials 29|. An energy 
cutoff of 40 Ry was used for the plane-wave basis set and 
integrations over BZ were performed using samples of 
1x8x8 k-points for monolayers and bilayers and 4x8x8 
k-points for bulk structures. 

For monolayer and bilayer models, we used periodic 
boundary conditions along the three dimensions, with 
vacuum regions of 8 and 9 A, respectively, between ad¬ 
jacent images in direction perpendicular to the layers. 
Convergence tests with greater vacuum spacing, guaran¬ 
tee that this size is enough to avoid spurious interaction 
between neighboring images. 

The optical conductivity was calculated directly from 
the joint density of states ie. taking into account only 
direct excitations. The real part of the dielectric function 
is then calculated using the Kramers-Kronig relationship. 
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FIG. 1. (Color online) Optimized structures of monolayers 
of group-IV monochalcogenides with phosphorene-like struc¬ 
ture. (a) Side view of the x-z plane for the four compounds 
and for phosphorene. (b) Side view of the y-z plane of SnS 
and phosphorene. (c) Top view of the structures, with the 
lattice vectors a and b along the x and j/-directions. (d) The 
respective BZ and the high symmetry points T, X, T and Y. 


III. RESULTS 
A. Crystal Structure 

The bulk a phase has an orthorhombic structure with 
eight atoms per primitive unit cell, four of each species. 
The primitive unit cell contains two puckered layers, 
stacked on top of each other. The bilayer is obtained 
by increasing the lattice supercell vector perpendicular to 
the plane of the layers. The monolayer has four atoms per 
unit cell, as in Fig.jljc). Each atomic species is covalently 
bonded to three neighbors of the other atomic species, 
forming zig-zag rows of alternating elements. Thus, there 
is in each atom a lone pair pushing its three bonds to¬ 
wards a tetrahedral coordination, just like in black phos¬ 
phorus, resulting in its characteristic waved structure. 

We adopt the axes system used by previous works [Tsi 
HO, mi, where layers are chosen to sit on the x-y plane 
ie. perpendicular to the z-direction. This is the same 
system conventionally used for black phosphorus, which 
takes the y-axis to be parallel to the puckering direction. 
Each layer has three non-trivial symmetry operations. 


namely a vertical mirror plane parallel to the xz plane, a 
two-fold screw rotation along an axis parallel to y, and a 
glide reflexion on a plane parallel to the xy plane. The 
atomic positions, in units of the unit cell vectors a, b 
and c, are ±(x,|,z; ^ -I- — x). For bulk SnS, for ex¬ 

ample, our calculated fractional atomic coordinates are 
a;(Sn)=z(Sn)=0.12, a;(S)=0.48 and z(S)=0.85. The lat¬ 
tice parameters and fractional atomic positions do not 
deviate much from those in the parent black phosphorus 
structure. As can be seen in Fig. UKa), the most no¬ 
ticeable difference is that the height of the atoms (along 
z) is no longer constant, rather cation and anion have 
slightly different heights alternating along x. The lattice 
is also more compact along a;, but less compact along the 
^/-direction, to decrease the repulsion between the atoms 
of the same type aligned up along the ripples (Table H]) . 
Full details of the structure of the other compounds, in¬ 
cluding the calculated a;(M,X) and z(M,X), can be found 
in the Supplemental Material (SM). 

The calculated lattice parameters are in good agree¬ 
ment with experimental data for bulk SnS (a=4.33A, 
b=3.98A and c= 11.20A) [1, [s^ and bulk SnSe (a=4.44 
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Monolayer 

Bilayer 


Bulk 


a 

b 

a 

b 

a 

b 

c 

SnS 

4.24 

4.07 

4.28 

4.05 

4.35 

4.02 

11.37 

SnSe 

4.36 

4.30 

4.42 

4.25 

4.47 

4.22 

11.81 

GeS 

4.40 

3.68 

4.42 

3.67 

4.40 

3.68 

10.81 

GeSe 

4.26 

3.99 

4.31 

3.97 

4.45 

3.91 

11.31 

P 

4.60 

3.30 

4.57 

3.31 

4.57 

3.51 

11.69 


TABLE I. Optimized lattice vectors in A for a phase of SnS, 
SnSe, GeS, GeSe along with those of phosphorene (P). 


A, b=4.15 A, and c=11.50).[33 Our results are also in 
agreement with previous theoretical studies of SnS [TTI . 
[l4 . . It is interesting to note that the lattice param¬ 

eters show little variation amongst the four compounds, 
differing less than 7%. This is due to the similar elec¬ 
tronegativity of Se and S, and of Ge and Sn, and conse¬ 
quently to the similar bond strengths. Thus, the lattice 
parameter trend is mostly dominated by the ionic radius 
of the constituents, the most compact structure being 
that of GeS. The lattice parameters of the sulphides re¬ 
main nearly unchanged for different number of layers, 
while those of the selenides show variations of ± 2%. 


B. Electronic Properties of monolayer, bilayer and 
bulk models. 

A tunable bandgap energy within the visible range, is 
one of the most interesting properties of the group-IV 
monochalcogenides. The calculated values for the energy 
gaps {Eg) are summarized in Table Hill These were ob¬ 
tained by calculating the bandstructures using the HSE 
functional along the high-symmetry paths of the Bril- 
louin zone (BZ) (FiglT]-d). The HSE exchange-correlation 
functional opens the gap, comparing to the PBE band- 
structure, while the band dispersion remains nearly un¬ 
changed 

Some aspects of the bandstructure are common to all 
systems studied. The dispersion of the bands nearest to 
the gap is nearly the same along the T-X and T-Y direc¬ 
tions, despite the striking difference between the struc¬ 
ture along those two crystallographic directions. This 
is probably due to the Sn-5s (Ge-4s) character of those 
bands. Additionally, in most cases there are multiple va¬ 
lence band and conduction band valleys, and most of the 
compounds have indirect gap, except for monolayer SnSe, 
bulk GeS and monolayer and bilayer GeSe, as shown by 
our calculated bandstructures in Fig. [2j This is different 
from black phosphorus and phosphorene, which has a well 
defined direct or nearly-direct gap. Finally, as expected, 
Eg decreases with increasing number of layers. 

In the following, we examine in detail the bandstruc¬ 
ture of each one of the compounds. The energy gaps 
given were obtained with the HSE functional, unless oth- 




Monolayer 

Bilayer 

Bulk 

SnS 

GBM 

(0.00 ; 0.84) 

(0.00 ; 0.74) 

(0.00 ; 0.70) 

VBM 

(0.22 ; 0.00) 

(1.00 ; 0.00) 

(0.12 ; 0.00) 

SnSe 

GBM 

(0.20 ; 0.00) 

(0.00 ; 0.78) 

(0.00 ; 0.72) 


VBM 

(0.20 ; 0.00) 

(0.20 ; 0.00) 

(0.30 ; 0.00) 

GeS 

GBM 

(0.00 ; 0.82) 

(0.00 ; 0.76) 

(0.00 ; 0.00) 


VBM 

(0.26 ; 0.00) 

(1.00 ; 0.00) 

(0.40 ; 0.00) 

GeSe 

GBM 

(0.20 ; 0.00) 

(0.20 ; 0.00) 

(0.80 ; 0.00) 


VBM 

(0.20 ; 0.00) 

(0.20 ; 0.00) 

(0.30 ; 0.00) 


TABLE II. (Golor online) (A*,,,, ; Aky)'. Positions of the VBM 
and GBM along the P-X and P-Y lines in the BZ. The values 
are given in units of 27r/|a| and 27r/|b| for ka, and kj, directions, 
respectively. 


erwise stated, and without considering spin-orbit cou¬ 
pling. Last, we will consider the spin-orbit splitting. 

a. SnS has indirect gap independently of the num¬ 
ber of layers. The indirect bandgaps calculated with HSE 
are Eg = 1.96, 1.60 and 1.24 eV for monolayer, bilayer 
and bulk SnS, respectively. The results for bulk agree 
very well with the experimental measured gap energies 
compiled in Ref. [l^, which presents bandgap energies 
obtained from optical absorption measurements and ex¬ 
trapolated to zero temperature around 1.2-1.3 eV, agree¬ 
ing very well with our calculated values. Another pre¬ 
vious study [sslj found a bandgap of 1.11 eV for bulk 
using HSE06, but different from our study, they used 
the experimental lattice parameters as input. The ex¬ 
cellent agreement between the calculated HSE bandgap 
and the experimental energies is unexpected, since the 
gap measured in optical experiments differs from the con¬ 
duction gap by the exciton binding energy, which can 
be of the order of hundreds of meV in two-dimensional 
materials. [1] However previous theoretical studies us¬ 
ing the GW method [Sy, [s^, obtained Eg = 2.57, 1.57 
and 1.07-1.26 eV for monolayer, bilayer and bulk SnS, 
respectively. [Ill EBl- This variation in the GW gaps 
may be due to the difficulty in treating the screening in 
2D materials. [1^ The agreement is very good except for 
monolayer, and this difference may be due to several fac¬ 
tors. One of them is the presence of shallow core d bands 
in SnS, as discussed in Ref. El- 

We now turn to the details of the bandstructure. In 
all cases (monolayer, bilayer and bulk), the valence band 
maxima (VBM) and conduction band minima (CBM) are 
located along the T-X and T-Y lines. In the monolayer, 
there are other competing local CBM and VBM, very 
close in energy to the band edges. By considering this, 
in addition to the indirect band gap observed for the 
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FIG. 2. (Color online) Electronic band structures for monolayer, bilayer and bulk group-IV monochalcogenides calculated using 
the HSE hybrid functional. The VBM and CBM are highlighted by full circles. Dashed black arrows indicate possible direct 
transitions (Ti and T 2 ) to points very close in energy to the VBM and CBM. Triangles indicate the position of the CBM when 
spin-orbit coupling effects are considered. 


monolayer, two direct gaps higher in energy by 75 meV 
(represented by transition Ti in Fig. [2]), and 0.22 eV (T 2 ) 
are also identified along F-X and F-Y lines, respectively. 
For the bulk, we also observe a competing point along the 
F-Y line, defining a direct gap of 1.40 eV. The situation 
is similar for the other compounds studied. 


b. SnSe For monolayer SnSe, the direct gap of 
1.44 eV is calculated along the F-X line. A second max¬ 
ima at 0.16 eV above VBM is obtained in the F-Y direc¬ 
tion, defining an additional direct transition at 1.60 eV 
(Ti). The bilayer bandstructure shows some differences. 
The CBM is now in the F-Y direction, and the material 
is characterized by an indirect gap of 1.20 eV, with the 
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Monolayer 


Bilayer 



Bulk 


GGA 

QE 

GGA 

VASP 

HSE 

VASP 

GGA 

QE 

GGA 

VASP 

HSE 

VASP 

GGA 

QE 

GGA 

VASP 

HSE 

VASP 


SnS 

1.40 

1.38 

1.96 

1.14 

1.12 

1.60 

0.83 

0.82 

1.24 

1.20-1.37 [H 

SnSe 

1.01* 

0.96* 

1.44* 

0.79 

0.76 

1.20 

0.55 

0.54 

1.00 

0.898 0, 0.95 [al 

GeS 

1.69 

1.65 

2.32 

1.55 

1.55 

2.20 

1.24 

1.22 

1.81 

1.70-1.96 [H 

GeSe 

1.14* 

1.18* 

1.54* 

1.02* 

0.98* 

1.45* 

0.59 

0.57 

1.07 

1.14 [13] 

P 

0.90 

- 

1.66 

0.55 

- 

1.30 

0.07 

- 

0.39 

0.33 [ilHil 


TABLE III. Gap energies {Eg) for monolayer, bilayer and bulk monochalcogeirides and phospliorene from GGA and hybrid 
functional calculations. We compare results from Quantum Espresso (QE) and VASP codes for the GGA approach. The star 
(*) indicates direct band gaps. Experimental values from previous works are also shown. All values are given in eV. 


VBM located in the F-X direction. 

For the bulk model, the CBM is also along the F-Y di¬ 
rection, and an indirect gap of 1.00 eV is defined with the 
VBM in the F-X direction. In this case, the second point 
nearest in energy to the VBM, lower in energy by 0.17 eV, 
is located along F-Y, and defines a direct gap of 1.17 eV 
(Ti in Fig. [2|). Both indirect and direct calculated gaps 
agree very well with the experimentally estimated values 
of 0.95 and 1.15 eV from Ref. [3^ . 



Monolayer 

Bulk 


Ti 

T2 

Ti 

SnS 

2.03 

2.18 

1.40 

SnSe 

1.60 

- 

1.17 

GeS 

2.62 

- 

1.83 

GeSe 

- 

- 

- 


TABLE IV. Smallest direct gaps obtained from the HSE cal¬ 
culations. The respective transitions are indicated by black 
arrows in Fig. [2] The values are given in eV. 

c. GeS Monolayer GeS has an indirect band gap of 
2.32 eV. This is defined by the CMB and VBM along 
the F-Y and the F-X lines. The lowest energy direct 
transition is at the F point (Ti in Fig. [2|). This direct 
gap is only 0.3 eV higher in energy than the indirect 
gap. In bilayer GeS, the gap is also indirect and 0.12 eV 
lower in energy than that of monolayer. The CBM in this 
case is still along F-Y, but the VBM is at the F point. A 
second VB maximum at the X point is almost degenerate 
in energy with the VBM at F. 

Bulk GeS has an indirect gap of 1.81 eV. The direct gap 
at the F point is only 20 meV higher in energy than the 
indirect one. Our results are in excellent agreement with 
the bulk GeS gap energies extrapolated to zero tempera¬ 
ture given in Ref. [l^, which range from 1.70 to 1.96 eV 
for different experiments. 

d. GeSe For GeSe monolayer, our calculations pro¬ 
duce a direct gap £'g=1.54 eV along the F-X line. For 
the bilayer structure, a direct gap of 1.45 eV is found in 


the F-X direction, near the X point. 

Similar to mono- and bilayer GeSe, the electronic 
structure for the bulk also shows the VBM and GBM in 
the F-X direction. However, the calculated gap of 1.07 eV 
is indirect in this case, as the GMB is located closer to 
the F point. In this case our calculations are also in very 
good agreement with experimental results, where optical 
measurements have indicated bulk GeSe as an indirect 
band gap semiconductor with £’g=1.14 eV [40j |. 


1. Spin-orbit splitting 

In the case of monolayer, the spin-orbit coupling gives 
rise to a splitting of the bands. This is different from 
monolayer phosphorene, where the inversion symmetry, 
together with time reversal symmetry, requires that the 
bands for the two spins are degenerate. No spin-orbit 
splitting is expected in bulk and even-numbered layer 
group-IV monochalcogeirides for the same reason. 



FIG. 3. (Golor online) Calculated SO-splittings of the CB and 
VB for the valleys along F-Y in meV. The schematic pictnre 
shows in continuous black lines the bands without spin-orbit 
effects, which are taken into account in a fully relativistic cal¬ 
culation. The lifted conduction and valence bands are shown 
in dashed lines. The lowering in Eg is given by A|V 

In order to quantify the effect of SO coupling, a fully 
relativistic calculation based on GGA was performed for 
the monolayer of all compounds. Fig. |4] shows the cal¬ 
culated electronic bands for monolayer SnSe with and 
without spin-orbit coupling. The results for the other 
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materials, given in SM, are very similar. So, for sake of 
simplicity, we discuss only the SnSe case. 

The spin degeneracy of the bands is lifted in all BZ 
except along F-X, the direction along which the C 2 rota¬ 
tional symmetry and the xz mirror symmetry are pre¬ 
served. Otherwise the SO-coupling changes little the 
shape of the bands, except for band crossings avoided in 
the relativistic result. However, the absolute minimum 
of the conduction band is along F-Y for all systems. For 
SnSe, the conduction band splits by 52 meV near the Y 
point (at the F-Y line). This defines a new CBM, since 
the other maxima at the F-X direction shifts by only 
38 meV. The same occurs for GeSe, for which a split of 
50 meV for the valley near Y is calculated. The VBM 
remains along F-X in all cases. The calculated spin orbit 
splittings are large for the conduction bands of all the 
phosphorene analogues, with the largest calculated value 
being for the SnS CBM (86 meV), as can be seen in the 
table presented in Fig. [3] This exceeds the spin-orbit 
splittings for the conduction bands of some of the most 
used transition metal dichalcogenides, eg. M 0 S 2 , MoSe 2 , 
WS 2 and WSe 2 , which have spin-orbit splittin gs b etween 
3 and 30 meV at the conduction band valleys, How¬ 
ever, the SO split on the valence band valleys is always 
smaller than for TMD, which have splittings of the order 
of 0.15-0.50 eV.[ii,|i| 



FIG. 4. (Color online) Electronic band structures for SnSe 
monolayer with (green and blue dashed lines) and without 
(continuous black lines) spin orbit coupling effect. 


C. Optical Properties 

We now analyze the optical properties of this class 
of materials, discussing the influence of anisotropy and 
the consequences of lower dimensionality. The complex 
dielectric function e(w) of the bulk materials has been 
measured by electron energy loss spectroscopy (EELS), 
optical transmission and reflectance measurements. [471- 
l50l | It can also be obtained from the optical conductivity 
tT(a;), which is related to e by: 


e(a;) = 1-I-cr(w), ct(w) = —fweo [e(w) — 1] • (1) 

weo 

where ui is the frequency of the incoming electromagnetic 
wave and cq the vacuum permittivity. 

The imaginary part of the dielectric tensor ei{u})a^i 3 for 
the bulk system can be obtained from the first-principles 
bandstructure using 


ij • 


M, 




{Ef{k) - E,{k)Y 


( 2 ) 


xd{Ef{'k) — Ei(k) — 


where £’i(k) represents the Kohn-Sham eigenvalues, Mq,^ 
represents the squared momentum matrix elements, a 
and /3 are the crystal directions, and subscripts i and / 
correspond to initial and final states, respectively. V and 
m are the cell volume and the electron mass. Note that 
this includes only direct transitions. 

From Eq. the real part of e{uj) can be obtained via 
the Kramers-Kroning relations: 


2 

er(w)ct/3 = 1-1-/ 

Jo 


uj'ei{uj')ap 

+ UJ^ 


du)' 


(3) 


Since the integral over the BZ needs a very large num¬ 
ber of points to converge, we use the GGA functional for 
the calculation of Ei{k),Ef{k), but subsequently apply 
a rigid shift to correct the bandgap to the value obtained 
from the HSE calculation. 

We start by validating the calculation by comparing 
tT(a;) obtained for bulk with the respective values ex¬ 
tracted from experiment. For all the four monochalco- 
genides, the agreement is excellent, as shown in Fig. [5l 
which compares the calculated real part of the conductiv¬ 
ity with experimental results from previous works. 
Both the threshold energy and the bandwidth of the at 
spectrum, which extends to about 20 eV in all cases, are 
well reproduced by the calculation. The low frequency 
value of the real part, Ur, is also in reasonable agreement 
with experiment. 

A striking feature is that the dielectric function is 
nearly isotropic i.e., the three diagonal components cr^x, 
(jyy and (Jzz have nearly the same magnitude and spec¬ 
tral dependence. This is a direct consequence of the sim¬ 
ilar band dispersion along the three respective directions 
of the BZ, which in turn originates in the structure. As 
shown in Fig. [U the Pcmn structure can be regarded as a 
distorted cubic NaCl structure, and the dispersion of the 
bands near the gap, mostly with s character, is little af¬ 
fected by the symmetry breaking, specially in plane. The 
nearly isotropic behavior of e is in agreement with previ¬ 
ous calculations, EM HH but at variance with experimen¬ 
tal data for GeS[47j. This may be due to the difficulty 
in performing measurements with the electric field polar¬ 
ization perpendicular to the layers. The nearly-isotropic 
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E(eV) 


FIG. 5. (Color online) Real part of the optical conductiv¬ 
ity ar((jj) for monolayer (full red lines), bilayer (full orange 
lines) and bulk (full blue lines) monochalcogenides. For com¬ 
parison, experimental results are also included from previous 
works. As the experimental works report ei(^, we use rela- 
tions[T]to obtain ar(cu). Ref(a): [i^ : Ref(d): j^; Ref(e): [4^ : 
Ref(f): [13. 


behaviour of the optical response is in variance with black 
phosphorus, which has significantly different absorption 
thresholds for the two in-plane directions. 

Plasmon energies of these materials can also be esti¬ 
mated from the first-principles dielectric function, when 
i?e[e] = 0. A comparison of the calculated real part of 
the dielectric function 6 ^( 0 ;) with experimental results 
from previous works is presented in Fig. [B] The plas¬ 
mon energies are calculated for the intrinsic (insulating) 
material ie. inter-band plasmons are considered. The 
calculated values are in the range 15-20 eV, in excellent 
agreement with experiment, [13 except for the sulphides, 
(Table |V| for which the plasmon frequency is slightly 
overestimated. We also note that although experimen¬ 
tally the plasmon frequency is the same for the three 
dimensions, there is a slightly variation in the theoretical 
results up to the uncertainty of the calculation. Even so. 


the experimental results are very well reproduced. 

We now compare the optical response of the mono- 
layer and bilayer material with those of bulk. For 
two-dimensional materials, the dielectric constant is not 
well defined, depending on the interlayer distance L 

asUllEllll 


e = 1-b 


47rx2D 

L 


(4) 


The 2D polarizability X 2 D is constant, and can be ob¬ 
tained from the first-principles calculations, where e is 
calculated from Eq. [2] by integrating over the whole su¬ 
percell. The interlayer distance L in that case thus corre¬ 
sponds to the supercell length along the direction perpen¬ 
dicular to the layers. [3 Hence, direct comparison between 
the 2D and 3D systems is not possible unless L is defined. 
Here, in order to obtain effective monolayer and bilayer 
dielectric constants comparable with bulk, thus high¬ 
lighting the differences originating in the bandstructure 
and lower dimensionality, we take L to be the bulk inter¬ 
layer distance |a|/ 2 , thus defining Cefi = 1 -|- (e — 1 )^- 
Comparing the effective optical conductivity for mono- 
layer, bilayer and bulk, we notice that the spectral de¬ 
pendence and magnitude are very similar, however, as the 
number of layers decreases, the peaks get sharper, due to 
the divergence of the joint density of states for saddle 
points of Ef(k) — Ei(k).[53 The real part of e is also 
nearly unchanged for effective interlayer spacing equal to 
that of bulk. More details of the points of interest in the 
spectra are given in Table IVl 


IV. CONCLUSIONS 

We have performed a systematic study of the electronic 
and optical properties of the family of phosphorene ana¬ 
logues SnS, SnSe, GeS and GeSe, and explored the conse¬ 
quences of lower dimensionality and symmetry breaking. 

One of the most interesting facets of these materials 
is that, from bulk down to monolayer, they cover a wide 
range of bandgap energies, from ^ 1.0 to 2.3 eV according 
to hybrid functional calculations. The bandgap increases 
as the number of layers is reduced, thus making it is 
possible to extend the absorption edge up to the green 
region of the spectrum. 

In parallel, this family of materials has the advantage 
of showing little variation in lattice parameters. This 
makes alloying a very promising direction for tuning the 
optical and electronic character of the materials. 

Further, the smallest lattice mismatch between SnS 
and GeSe compounds, for both lattice parameters, can 
be an indication of these materials as good candidates 
for formation of hybrid structures, [s^ Phosphorene itself 
has a small lattice mismatch to some of these materials, 
with which it can be combined. 

Additionally, the inversion symmetry breaking in 
monolayer allows for spin-orbit splitting of the conduc¬ 
tion and valence band valleys along the F — F direction. 
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OJp 




'-r 


LiJp 


tUp-Exp. 


X 

y 

X, y 

X 

y 

X 

y 

z 

x,y,z 

X 

y 

z 

x,y,z 

SnS 

9.9 

10.0 

0.69 

19.3 

19.2 

10.30 

10.94 

9.90 

0.72 

18.1 

17.8 

18.0 

16.22 

SnSe 

12.5 

12.8 

0.77 

14.8 

14.8 

10.90 

11.50 

9.90 

0.78 

15.7 

15.7 

15.9 

15.44 

GeS 

8.7 

8.6 

0.63 

20.4 

20.4 

9.20 

9.45 

9.10 

0.64 

20.2 

20.1 

20.2 

18.21 

GeSe 

13.8 

14.7 

0.70 

17.2 

16.5 

10.60 

10.96 

9.80 

0.72 

17.7 

17.2 

17.8 

17.30 


TABLE V. Calculated values of the real part of the dielectric constant at tu —>■ 0 (e°) and at high energy limits The 

plasmon frequencies (ojp) for x and y-directions in monolayer and x, y and 2 -directions in bulk are also presented. For both 
monolayer and bulk, £“ assume the same values for all directions, while the static dielectric constant ejl and the plasmon 
frequency ujp varies slightly for different directions. Experimental values of ujp for the bulk of all compounds from Ref. [4^ are 
also included. All values are given in eV. 



E(eV) 



E (eV) E (eV) E (eV) 


FIG. 6. (Color online) Real part of dielectric constants er(w) 
for monolayer (full red lines), bilayer (full orange lines) and 
bulk (full blue lines) monochalcogenides. Experimental re¬ 
sults from previous works are also included. Ref(a); [i^ : 
Ref(b): [13]; Ref(c): [sl; Ref(d): [13]. 


This effect is absent in bulk group-IV monochalcogenides 
and in phosphorene. The spin-orbit splitting can be as 
large as 86 meV for the conduction band minimum of 
SnS. The large spin-orbit splitting is likely due to the s 
character of the VBM and CBM. 

The properties described here open the possibility of 
using group-IV monochalcogenides for optoelectronics 
and spintronics. Despite, to our knowledge, there being 
still no experimental reports of isolation of monolayers of 
the monochalcogenides considered here, few-layers isola¬ 
tion by exfoliation has already been achieved pus, 
and it is just a matter of time until control of the layer 
number becomes possible. 

Finally, this illustrates that phosphorene is the parent 
structure of a whole class of orthorhombic 2D materials, 
whose potential is to be revealed once they are isolated 
in monolayer or few-layer form. 
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